Second generation of vortex-antivortex states in mesoscopic superconductors: 

stabilization by artificial pinning 
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Antagonistic symmetries of superconducting polygons and their induced multi- vortex states in a 
homogeneous magnetic field may lead to appearance of antivortices in the vicinity of the supercon- 
ducting/normal state boundary (where mesoscopic confinement is particularly strong). Resulting 
vortex-antivortex (V-Av) molecules match the sample symmetry, but are extremely sensitive to 
defects and fluctuations and remain undetected experimentally. Here we show that V-Av states 
can re-appear deep in the superconducting state due to an array of perforations in a polygonal set- 
ting, surrounding a central hole. Such states are no longer caused by the symmetry of the sample 
but rather by pinning itself, which prevents the vortex-antivortex annihilation. As a result, even 
micron-size, clearly spaced V-Av molecules can be stabilized in large mesoscopic samples. 

PACS numbers: 74.78.Na, 74.25.Dw, 74.25.Qt 
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Vortices in small superconducting elements have been 
the focus of scientific research ever since it was shown 
that a strongly confined superconducting condensate 
shows enhanced critical properties 0,0. Usually indi- 
vidual vortices pierce the sample in the direction of ap- 
plied magnetic field, minimizing the energy of the stray 
magnetic field, as well as their mutual interaction. Such 
formed vortex clusters (called 'multi- vortex' tend to 
mimic the symmetry of the sample, due to the repulsion 
with the surrounding screening (Meissner) currents. In 
samples of the order of the coherence length £ and/or 
penetration depth A, the edge currents may even com- 
press the vortex lines into a single bundle, often referred 
to as a 'giant- vortex' However, in addition to these 

simple arguments, the effects of lateral confinement in 
mesoscopic elements can have more elaborate manifesta- 
tions. One such is the appearance of vortex-antivortex 
molecules, of exact the same symmetry as the host sam- 
ple. 

Almost a decade since the original prediction of Chib- 
otaru et al. [5| the symmetry-induced antivortex remains 
one of the most exciting phenomena in vortex matter in 
submicron superconductors. The stability of an antivor- 
tex in opposite magnetic field is still puzzling, and not yet 
verified experimentally. Latter is not surprising, taking 
into account that predicted vortex-antivortex molecules 
are very small in size (~ £, somewhat larger in type-I 
samples [(|) and extremely sensitive to defects Q- Actu- 
ally, in most instances, they are simply unstable @, 0|. 

In last several years, two methods were proposed to 
improve the stability and the observability of the vortex- 
antivortex states in mesoscopic polygons. In one, the 
magnetic field profile was altered by a magnetic dot 
placed on top of the superconductor [nj. The bipolar 
field of the added magnet favors the antivortex under- 
neath, and repels vortices further apart. However, this 
method: (i) interferes with the key concept of an antivor- 
tex in opposite, unipolar field, (ii) makes the structure 



more complicated, three-dimensional, and inaccessible 
for scanning probe techniques, and (iii) poses difficulties 
for any magnetic (magnetic-force, Hall-probe) measure- 
ments. Instead, we recently proposed the structural engi- 
neering of the sample itself, by strategically placed holes 
in the sample, mimicking the sample symmetry and the 
expected symmetry of the vortex-antivortex state ll| . In 
such realization each hole hosts (pins) one vortex, where 
the pinning force is effectively stronger than the vortex- 
antivortex attraction. As a consequence, vortices remain 
in the holes even when placed at further distances from 
the antivortex, and the whole molecule can be made sig- 
nificantly larger. 

In this paper, we realized that not only vortices can 
be pinned in the latter concept, but an antivortex as 
well. We therefore introduced an additional, central hole 
in the sample, in which the antivortex can reside. Intu- 
itively, this facilitates the quantization of negative stray 
flux between the vortices in a cluster, and leads to a more 




FIG. 1: A schematic drawing of a polygonal mesoscopic super- 
conductor (side a, thickness d) with corresponding arrange- 
ment of perforations (spacing In, characteristic size of each 
hole ah and of the central hole a c ) in a perpendicular mag- 
netic field H. 
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stable antivortex. Vortices and antivortices in holes are 
just quantized magnetic field, screening prevents their 
annihilation, and the fact that they have no core reduces 
their attractive interaction. As will be shown, this in- 
deed leads to an enhanced stabilization of the vortex- 
antivortex molecule, even deep into the superconducting 
state, but with significantly altered overall properties. 
Although this generally applies to any polygonal setting 
of the holes and various shapes of the sample 
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what follows we show the proof of principle for a square 
superconducting geometry, as shown in Fig. [TJ 

For the purpose of this study, we us the mean-field, 
Ginzburg- Landau (GL) theory, and solve numerically 
the set of GL equations 



(-»V - AYi> — (1 — T - ivr)?A = 0, 
-k 2 AA 



j s = $ V UV - A) tp* 



(1) 
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already established as a prime tool for theoretical de- 
scription of mesoscopic superconductors Q. All vari- 
ables in Eqs. (|H2p are dimensionless and expressed in 
temperature independent units (assuming ~ (1 — T)^ 1 / 2 
dependence of characteristic length scales in the system, 
coherence length £, and penetration depth A), n = A/£ 
is the so-called GL parameter, and describes the ability 
of the sample to screen the applied magnetic field. Tem- 
perature is scaled to T c , magnetic field to H c2 (0), vector 
potential to if C 2(0)£(0) and the order parameter ip to 
ip(0) = y/a(0)/{3. The Gibbs free energy density T , in 
units of JTT=0, H=0)= H c (0) 2 /8ir is calculated from 
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HV-A> -(l-T-^\ 2 M 2 (3) 
+k 2 (V x A- H) 2 ]dV, (4) 



for each of the stable superconducting states found dur- 
ing the simulation. The full simulation region, used for 
the calculation of the magnetic response of the sample, 
was a square with 128 grid points in each direction. The 
number of grid points inside the sample was kept at 
64 x 64. 

In Fig. [5] we show the free energy of vortex states 
with vorticities from L = 0toL = 6asa function of 
the applied magnetic flux $ = a 2 H . The simulation is 
performed by sweeping the magnetic field up and down, 
and then back-tracking of all found vortex states so that 
their complete stability range is obtained. In addition, we 
performed the calculation starting from different initial 
conditions, some of which included reasonable guesses of 
possible vortex configurations. Fig. [3 summarizes all the 
found stable states, but we may not rule out the possibil- 
ity of some more complex, higher energy, vortex states. 
For this simulation a square of size 20£(0) x 20£(0) was 
used with holes of size 3.5£(0) x 3.5£(0), four of which are 
horizontally and vertically displaced by 5£(0) from the 
center of the sample. Temperature was fixed at 0.3T C 



(corresponding to Al samples below 400 mK [13j) and 
an extreme type-II behavior was assumed, typical for 
thin samples (k — oo). The combinatorial number of 
possible vortex configurations for given vorticity is quite 
high, and many of them are indeed found stable due to 
both the large size of the sample and that of the cen- 
tral hole. Starting from zero field, the Meissner state, 
L = 0, is the lowest energy state in the flux band of 
A$ = 1.85<!>o, which is larger than a flux quantum, but 
significantly smaller than in non-perforated mesoscopic 
samples [9]. The reason is that the first vortex pene- 
trates easier, and is more stable in the present system 
due to the central hole (central L = 1 is in the ground 
state in A$ = 1.92$o). States with an off-center vortex 
have mostly higher energy, smaller stability range, and 
are fourfold degenerate (which makes them interesting 
for logic applications [IH). L = 2 shows the broken sym- 
metry in the ground state, as two vortices reside in the 
central and corner hole, but their ground-state flux-band 
is only A$ = O.lT^o- Similarly, the broken symmetry 
L = 3 state with all vortices in the corners is the ground 
state in a short interval of A$ = 0.11<i>o. Just like in 
the superconducting square with 2x2 holes [l5|], L = 4 
vortices reside in the corners and show enhanced stabil- 
ity in the ground state (A<i> = 1.95<I>o). However, the 
most dominant state in Fig. [2] is the L — 5 one, where 
all holes are occupied by a single vortex (A$ = 3.63$o)- 
Therefore, while for the square with an 2 x 2 holes high 
symmetry states as L = 2 and L = 4 were pronouncedly 
low-energy states, we notice that in the present system 
the states which have a vortex in the central hole, and 
are fourfold symmetric dominate the ground state, i.e. 
L = 1 and L = 5. Note that the L — 4 state is still very 
comfortable with the fourfold symmetry although it does 
not have a central vortex, but e.g. L = 3 with a central 
vortex is not in the ground state because it does not obey 
the symmetry of the sample. Because of their very low 
energy, L = 1,4,5 states overshadow other states like the 
L = 2 and L = 3. The L — 3 state has several allotropes, 
i.e. it can have vortices across the sample diagonal, one 
vortex in the center and two on the side, or all three out 
of the center. Only one configuration can satisfy the four- 
fold symmetry of the sample and that is the L = 4 — 1 
state, with four vortices in corner holes and an antivor- 
tex in the center. This vortex-antivortex (V-Av) state is 
indeed found stable in our system, with relatively high 
energy, but still lower than most higher vorticity states 
(which suggests that it can be experimentally realized in 
decreasing magnetic field). 

The physical origin of this V-Av state is however differ- 
ent from the V-Av state found in Refs. 0, @, EH where 
the V-Av molecules appear due to the imposed symme- 
try of the sample. In Ref. [llj, the added holes only 
reinforce the symmetry argument and do not cause per 
se the V-Av state. In the rest of the paper, we will refer 
to this state as the symmetry induced V-Av in contrast 
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to the newly found V-Av state which is fully stabilized 
by pinning of all vortices and the antivortex. In Fig. [3] 
we show the $ — T phase diagram for both versions of 
the L = 3 vortex-antivortex state (for taken parameters 
a = 20£(0), a h = 3£(0), l h = 4£(0), a c = 1.5£(0), see 
Fig. Q]). We observed two fully independent islands in 
the $ — T phase space corresponding to the two mani- 
festations of the V-Av state. The symmetry-induced one 
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FIG. 3: The superconducting state in the $— T diagram (solid 
line shows the superconducting/normal state boundary) com- 
prising the areas of stability of the L = 3 vortex-antivortex 
molecules, both pinning- and symmetry-induced ones. Color 
coding shows the energy of those vortex states. 



is situated in the high temperature regime, where con- 
finement imposed by the sample boundaries is effectively 
very large. It should be noted here that the AT sta- 
bility range for the symmetry-induced V-Av state is on 
average just 1% larger than in the case of 2 x 2 holes 
(since the order parameter inside the molecule is in 
any case severely suppressed). At lower temperatures, 
the GL equations are strongly non-linear and symmetry 
arguments cannot account for the nucleation of the V- 
Av state. Nevertheless, a second generation V-Av state 
does stabilize, thanks to the large spacing between the 
holes (preventing vortex-antivortex annihilation) and the 
large size of the sample (diminishing the effect of encir- 
cling screening currents). The pinning-stabilized V-Av 
state: 1) occupies far larger $ — T area compared to 
the symmetry induced one, a wide region of ks 10$> and 
« 0.6T C , and 2) it is found stable even in negative applied 
field, where central antivortex is a natural state and vor- 
tices are subjected to an increasing expulsion pressure. In 
principle, several vortex states are found stable at nega- 
tive fields in Fig. [21 which is a manifestation of the flux- 
trapping effect [16]. However, out of all L = 3 states, 
L — 4 — 1 shows maximal resilience to negative flux. 

To emphasize again, the size of the sample is playing 
a crucial role for the V-Av state. While the symmetry- 
induced V-Av exists only for small samples [compared to 
the coherence length £(T)] the pinning-stabilized V-Av 
requires a large sample. This is illustrated in Fig. HJa) 
for a square with five holes. When the sample is larger, 
the flux stability region of the pinning stabilized V-Av 
state becomes larger as well. This we can attribute to 
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the weakened influence of the Meissner current on the 
inside vortices because of the larger distance. One also 
notices the extremely large flux interval (from — 7.5$o to 
17<&o) for stability of the L = 3 pinning-stabilized V-Av 
state. However, penetration field for new vortices must 
decrease for larger samples (see e.g. 17] for the case 
of large BiSCCO disks). It is already known that mag- 
netic field higher than H c i is needed for penetration of 
vortices into mesoscopic superconductors, but this factor 
decreases to unity in bulk systems. Nevertheless, lat- 
ter factor decreases with size of the sample much slower 
than the square power increase of the area of the sam- 
ple. As a consequence, the threshold flux for penetration 
of new vortices in our sample increases as a function of 
a/£(0). We also observed the changing curvature of the 
latter dependence (at a/£(0) « 25); from the calcula- 
tion of threshold magnetic field, we found that samples 
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FIG. 4: The stability diagram vs. applied field (i.e. magnetic 
flux) for the pinning-induced vortex-antivortex state, (a) as a 
function of the sample size for a fixed temperature, and (b) 
as a function of the size of the central hole. 
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FIG. 5: The energy of V-Av molecules as a function of temper- 
ature, for different values of the Ginzburg-Landau parameter 
k. In numerical simulations, an abrupt cooling enables direct 
transition between two kinds of V-Av states. 



with size of 10-25£(0) exhibit almost identical demagne- 
tization effects i.e. similar vortex penetration field which 
then gradually decreases for sizes above 25£(0). 

The size of the central hole is equally important as the 
sample size, as it relates to the ability of V-Av pair to 
annihilate. Consequently, this second generation V-Av 
state cannot be stabilized in the absence of the central 
hole [see Fig. 0Jb)]. The effect of decreasing the cen- 
tral hole size is also indicated in Fig. 0Jb) - the upper 
boundary of the flux stability interval is decreasing, due 
to facilitated V-Av annihilation (at high fields, vortices 
are pushed towards the center of the sample), while the 
lower boundary is in fact not influenced by the size of 
the central hole. Note that lower boundary lies in the 
negative field region, where vortices have expulsion ten- 
dency, and pinning of any strength stabilizes the central 
antivortex with ease. 

The latter stability of the pinned V-Av state in the neg- 
ative field forms a base for the realization of this state 
in experiment. Namely, abrupt increase of the magnetic 
field from a low negative value (that stabilizes one an- 
tivortex in the sample) to a large positive one (allowing 
for penetration of multiple vortices in outer holes) may 
result in the desired V-Av state at low temperatures. 
Alternatively, following the results of Fig. [3l one can 
think of an abrupt temperature decrease. We performed 
a corresponding simulation, where we started from the 
symmetry-induced V-Av state (which is the ground state 
for L = 3 at temperatures near T c ) and then applied 
a steep temperature decrease. Without taking into ac- 
count the experimentally relevant temperature gradient 
over time, the simulations were always able to land into 
the pinning-stabilized V-Av state at low temperatures. 
This is illustrated in Fig. [5] where also the influence of 
k is depicted. Stronger screening of the magnetic field 
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FIG. 6: Calculated magnetic field profiles in the sample: (a,b) 
for the pinning-stabilized V-Av state and two values of k, and 
(c) the corresponding profile for the symmetry-induced V-Av 
state. 



into the holes (i.e. lower k) enhances the pinning, and 
consequently favors the V-Av state of second generation, 
in contrast to the symmetry-induced V-Av which is dis- 
favored by low values of k due to the effective attraction 
of vortices into a single bundle (see Ref . [ll[ ) . 

Although quite difficult to realize (never being the 
ground-state of the system) , the second generation V-Av 
state is much more suitable for direct experimental vi- 
sualization than its symmetry-induced ancestor. Firstly, 
samples can be made larger, and spacing between vortices 
and the antivortex are sufficient for their identification. 
Second, the Cooper-pair density is significantly higher 
inside the V-Av molecule, since temperature may be far 
below T c . Third, the magnetic field profile in the sam- 
ple is sufficiently inhomogeneous, with pronounced am- 
plitudes, benefiting from the size of the V-Av state and 
its stability for low values of k. This is shown in Fig. [5] 
where the magnetic field profile in and nearby the sample 
is calculated for two values of k for the second genera- 
tion V-Av and for one value for the symmetry-induced 
V-Av (given in real units). For taken parameters, the 
magnetic field contrast inside the symmetry-induced V- 
Av molecule is of the order of 0.01 Oe, while that of the 
second generation V-Av can be above 10 Oe, thus three 
orders of magnitude larger. 

Conclusion Using a square superconducting geome- 
try with five holes, we found a second generation of the 
L = 3 vortex-antivortex (V-Av) state in a uniform field. 
This state is purely induced by pinning, and is NOT 
caused by the symmetry of the sample, contrary to the 
previously found V-Av molecules. This novel state is en- 
ergetically favored in rather large mesoscopic samples, 
but it is never the ground state of the system. Con- 



sequently, more elaborate techniques are needed for its 
stabilization in experiment, such as an abrupt increase 
of field (from negative to positive value), or an abrupt 
cooling of the symmetry-induced V-Av state. However, 
the very large $ — T stability range of the second gen- 
eration V-Av state, the very comprehensive size of the 
V-Av molecule, and the large variation in amplitudes of 
both superconducting order parameter and stray mag- 
netic field inside the molecule (further enhanced for lower 
k and temperature), all lead to facilitated experimental 
observation. 
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